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Abstract. The dependence of the individual and of correlated skeletal variables of Carduelinae and Fringillinae upon bo¬ 
dy mass (=size) was tested by means of multiple regression and single linkage R-cluster analyses. The most obvious re¬ 
sult is the significantly high degree of body mass dependences of all individual variables within the three functional com¬ 
plexes feeding, hind limb locomotion and flying. Elimination of body mass as a measure of size decoupled several units 
of correlated variables. Nevertheless the general robustness of units was strong but the removal of “size” had an impact 
on the extent of morphological integration, measurable as agglomeration level of variable units. Body mass (size) in our 
sample thus influenced the quantity rather than the quality of morphological integration. The key feature in fringilline/car- 
dueline skeletal morphology - coevolution of skull and hind limb traits - thus appears to be relatively unconstrained by 
size. 

Keywords. Size, ecomorphology, functional complexes, biological roles, skeletal morphology, developmental gene¬ 
tics. 


1. INTRODUCTION 

1.1. The role of body size in avian morphology 

Body size has been identified as an important factor driv¬ 
ing the diversification and evolution of organisms. Olson 
& Miller (1958, 1999), e. g. quantified size as a strong 
integrative force of morphology, Fairbairn ( 1997) iden¬ 
tified its role in sexual dimorphism, Klingenberg & 
Spence (1997) analyzed its function in life-history evo¬ 
lution. Bonner (2004) proposed a third size dependent 
rule besides the already well known “weight-strength-rule” 
and “weight-surface-area-rule”. His “weight-complexity- 
rule” explains that increase and decrease in size of organ¬ 
isms are strongly related to increase and decrease in com¬ 
plexity in general, from differentiation of cell types to the 
social organization of organisms. 

These examples implement size as a morphological key 
feature in the organization of life, at all hierarchical lev¬ 
els from the individual to species (Barbosa et al. 2000; 
Bonner 2004; Chernoff & Magwene 1999; Peters 
1993). 

Size, as the physical magnitude of objects, can be quan¬ 
tified using either measures or weights. Linear length 
measures are used to seize both the two dimensional area 
as well as the three dimensional volume, whereas weight 
quantifies the proportionality to volume alone. Measures 


widely accepted to represent overall size of organisms in¬ 
clude also factors like the first principal component 
(Bookstein et al . 1985; Zelditch & Fink 1995; Leisler 
& Winkler 2006), landmarks or ratios. Bjorklund & 
MerilA (1993), e.g. used orthogonal rotation to remove 
size from species mean vectors, Leisler (1980) calculat¬ 
ed the cube root of a compound measurement of the body 
core skeleton (sternum length plus pelvis length plus cora¬ 
coid length) times (sternum width plus pelvis width) times 
height of crista stemi as measure for “body size”, follow¬ 
ing Hoerschelmann (1966). 

For the characterization of volume, body mass is the ad¬ 
equate measure. Several rules explaining the construction¬ 
al constraints on morphology are related to body mass: The 
“weight-strength-rule” already mentioned above, ex¬ 
plains why e.g longer arms, wings and legs have to be 
proportionally thicker than shorter ones, the “weight-sur¬ 
face-area-rule” describes physiological relations of e.g. gas 
diffusion in large versus small lung surfaces. These are 
geometric growth laws, that explain the change of propor¬ 
tions - so-called allometric relations - within a biologi¬ 
cal organism. The underlying process is differential growth 
during development brought about by differences in the 
growth rates of the various parts (Waddington 1966). 

In small bodied flying vertebrates, flight imposes identi¬ 
cal physical constraints on the flying organism that lead 
to a functional symmetry in body sizes (Maurer et al., 
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2004). The body plan of birds is primarily constrained by 
the weight-strength-rule, that regulates Hying ability. In 
extant species 22 kg is the limit for take off; the world¬ 
wide tiniest bird is the Bee Hummingbird, Mellisuga he- 
leiiae, found in Cuba and weighing from 1,6-1,9 g (Chai 
1999). The most successful clade in modem birds, the 
passerines (or perching birds, Passeriformes) have minia¬ 
turized to an average body mass minimum of 4,5 g in some 
tropical species (Black-faced Flycatcher Warbler, Abro- 
scopus schisticeps, Broad-billed Flycatcher Warbler, 
Ticket li a hodgsoni, Cuban Gnatcatcher, Polioptila lem- 
beyei; Dunning 1993 ). The heaviest and thus largest song 
bird, the Raven, Corvus corax, from Alaska weighs 1.240 
g. Perching birds successfully evolved and radiated with¬ 
in any terrestrial habitat -with exception of the polar re¬ 
gions where their occurrence is physiologically con¬ 
strained by the weight-surface-area-rule. 

1.2. Key features of the carduelid skeleton 

Our paper focuses on the impact of body mass on the ar¬ 
chitecture of the avian skeleton, taking the song bird sub¬ 
families Fringillinae and Carduelinae as a case study. Mor¬ 
phological diversity in skeletal elements in this phyloge¬ 
netic entity (clade, monophylum) is mostly confined to 
changes in overall size (Bjorklund 1991, 1994; van den 
Elzen & Nemeschkal 1986, 2006; Nemeschkal & van 
den Elzen 1990, van den Elzen & Nemeschkal 2006). 
Differentiations in shape are (despite their unquestioned 
ecological importance) only expressed as an insignificant 
portion of the total variation. Whereas three shape axes 
explained only 9,3 % of variance, the size axis account¬ 
ed for 85,5 % of variance (van den Elzen & Nemeschkal 
2006). The observed pattern has been identified as a gen¬ 
eral rule in passerine birds by Bjorklund (1994, 2006). 

. ..making species to be larger or smaller copies of each 
other, all moving up and down a common line of allo- 
metry”. Mapping of skeletal variables along the first prin¬ 
cipal components (Nemeschkal & van den Elzen 1990) 
revealed, that morphological diversification in Fringilli¬ 
nae and Carduelinae focuses on variables of the beak. 
Thus, morphological change in these subfamilies is main¬ 
ly reflected in feeding mechanisms. 

In terms of ecomorphology, a discipline taking the bio¬ 
logical role of morphological structures under consider¬ 
ation (Bock & von Wahlert 1965), elements of the avian 
skeleton belong to three different functional complexes: 
skeletal elements of the skull and beak are allied to the 
feeding complex, wing bones, pectoral elements and bones 
of the shoulder girdle represent the flying apparatus and 
skeletal elements of pelvis and legs the hindlimb locomo¬ 
tion complex. In the carduelid skeleton, nevertheless, units 
of covarying variables (UCCs) contain combinations oth¬ 


er than proposed by the functional complexes described 
above. These UCCs were detected to reflect units carry¬ 
ing out common biological roles. Besides a more or less 
undisturbed module of the “flying complex”, formed sole¬ 
ly by wing bones plus elements of the shoulder girdle, the 
variables of the feeding and hindlimb locomotion com¬ 
plexes were united in one single common cluster, and 
demonstrate the superimposure of feeding on hindlimb lo¬ 
comotion. 

In search of additional modern approaches in the study of 
morphological integration, Chernoff & Magwene (1999) 
propose a hierarchical framework for integrative hypothe¬ 
ses: Morphological variables are judged as integrated at 
the broadest, i. e. the uppermost inclusive level due to their 
covariation with size. At a less inclusive level variables 
covary due to developmental and/or functional associa¬ 
tions and at still lower levels due to anatomical (spatial) 
associations. Starting out from physics we may assume 
that changes of singular parts from a macroscopical phys¬ 
ical (e. g., mechanical) system will be ruled by various 
law-dependencies when the whole system changes size 
(e.g., mass). As a consequence, we test in the current study 
the dependence of individual and correlated skeletal vari¬ 
ables upon body mass taken as a measure of size. 

1.3. Presumptions 

Based upon our knowledge from preceding analyses (van 
den Elzen, Nemeschkal & Classen 1987, Nemeschkal 
& van den Elzen 1990, Nemeschkal, van den Elzen & 
Brieschke 1992, van den Elzen & Nemeschkal 2006) 
our hypotheses are that: 

1. a significant dependence of skeletal variables on body 
mass exists. It is expected to differ between variables 
and to be highest in beak elements; 

2. thus the three theoretical functional complexes -feed¬ 
ing. hindlimb locomotion, flying- differ in their size de¬ 
pendence; 

a. variables of the feeding apparatus - showing 
highest morphological differentiation between 
species - should exhibit the largest amount of 
size dependence, 

b. elements of the hindlimb locomotion complex 
and of the flying apparatus are thought to show 
lesser dependence, as of lesser interspecific 
variability. 

3. units of variables (UCCs) are also expected to be size 
dependent, elimination of “size” is thought to decouple 
variables within units; decoupling is expected to occur 
particularly within “mixed” variable units of beak and 
hindlimb osteometries. 
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2. MATERIAL AND METHODS 

In 313 specimens, representing 43 species (Table 1), 42 
equivalent skeletal variables (measurements of bone 
lengths, widths and depths) were taken, as described in 
detail inNEMESCHKAL (1999). Twenty measurements stem 
from the feeding apparatus, 13 represent the flying appa¬ 
ratus and 9 hindlimb locomotion (Table 2). 

Firstly, centroids of 43 species over 42 variables were built 
from the log scaled original data matrix consisting of 313 
specimens. This procedure was essential, because body 
mass data were available as species means only. Morpho¬ 
logical variation between species was then quantified by 
the total variance of species centroids. Body mass varia¬ 
tion in the actual fringilline-cardueline sample ranges from 
a maximum of 54 g in the Hawfinch Coccothraustes coc- 
cothranstes to a minimum of 9 g in neotropical siskins 
(Table 1). To test hypotheses about size impact, linear re¬ 
gression analysis was applied. Log transformed species 
means of body mass were taken as predictor variable and 
each of the 42 skeletal measurements as criterion variable. 
The 42 resulting residuals are used as the variables under 
study - the variables corrected for body mass (= variables 
of which size was partialized out; BCM). The coefficients 
of determination between body mass and original variables 
(Table 2) were tested for significance using random boot¬ 
strap (1000 replicates each; for computer programme 


package see Nemeschkal (1999)). They are figured as 
profiles of variance (Figs 1, 2). Single linkage R-cluster 
analyses were chosen to figure correlations of original da¬ 
ta (Figure 3- BO; based upon the variable intercorrelation 
matrix between species means for original data) and cor¬ 
relations data corrected for body mass (Fig. 3- BCM; vari¬ 
able intercorrelation matrix between species means for 
body mass adjusted data). Units of correlated variables 
were taken into account, when they agglomerated at the 
uppermost quartile (r 2 > 0.924) in the original dataset and 
reappeared in the size corrected dataset again. E-units are 
stable clusters with topological identity of variable posi¬ 
tions in both the R-cluster analyses of original data and 
the body mass corrected data set. S-units are less stable 
units built of identical variables but with topological dif¬ 
ferences between analyses of original and size corrected 
variables. 


3. RESULTS 

The most obvious findings are on the one hand the sig¬ 
nificantly high degree of body mass correlation of all in¬ 
dividual variables (Fig. 1) and on the other hand the ro¬ 
bustness of several covarying variable units found in both 
the original and the “residual” data set, adjusted for body 
mass (Fig. 2). 



Fig. 1 . Dependence of 42 skeletal characters on body mass as shown by coefficients of detennination. 












variances and residuals 
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BO-BCM 


Fig. 2. Comparison of variance profiles. Variance of original characters (BO) occur in the upper line and residuals for characters 
corrected for body mass (BCM) at the bottom. 


0.696 


0.772 


0.848 


0.924 



1.0 


0.594 


0.391 



0.188 


Figs 3. Single linkage R-cluster analysis of original data (BO; left) and of body mass corrected data (BCM; right). Black squares 
indicate measures of the pectoral girdle and forelimbs, black triangles mark skull measurements and open circles measures of the 
pelvic girdle and hindlimbs. Clusters are explained in the text. 
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3.1. Influence of body mass on individual variables 

All 42 variables are correlated with body mass at a high 
significance level (p< 0.001) as revealed by random boot¬ 
strap procedures. The degree of correlation with “size”, 
measured by coefficients of determination (COD ) between 
(log transformed) original data of species means and body 
mass, is highest in coracoid length (variable number 35; 
COD=0,91), post orbital length (12; 0,90) and equal at 
COD= 0,89 in tibiotarsus width (27), total skull length (16) 
and synsacrum length (21). The largest amount of body 
mass independence exhibit carpometacarpus (41) and ul¬ 
na lengths (39; COD= 0,34 and 0,47) as well as tar- 
sometatarsus length (No 28) and width (29; COD= 0,70, 
resp. 0,67) and proximal end width of the humerus (38; 
COD= 0,71). Skull variables in general reveal the high¬ 
est correlation with body mass, hindlimb and forelimb 
variables indicate a larger variation in their correlation 
with body mass. 

A comparison of variance profiles of original and size cor¬ 
rected skeletal variables (Fig. 2) confirms and strength¬ 
ens these findings: The relatively higher variances in the 
original variables appear smoothened in the size correct¬ 
ed residuals (BCM; bottom line in fig. 2), indicating that 
the largest amount of total variance between species is due 
to a “size factor”, measurable by body mass. This espe¬ 
cially holds true for skull variables (skeletal variables 
1-20) and less for hindlimb and forelimb measures (skele¬ 
tal variables 24-29, 37^42). 


3.2. Influence of body mass on united variables 

The single linkage cluster analysis of original data (Fig. 
3) extracted two main units fused at a high level, several 
disintegrated variables and one isolated variable set (El). 
The module BOC 1 is composed of the variables of the 
flying apparatus, mainly its “engine-supporting” part, in¬ 
cluding measures of sternum, coracoid and scapula. Wing 
measures that represent “flight performance” are highly 
disintegrated, only humerus length (37) and carpometacar¬ 
pus width (42) are included. The second main module 
BOC 2 combines skull and hindlimb measures. 

The two main modules BOC 1 and BOC 2 are composed 
of two different types of subsets: Robustness is observed 
in nine smaller units (E1-E9), where only two to three 
variables are coupled. Similar structure is maintained af¬ 
ter correction for size in four units (S1-S4), including three 
to 17 variables as well as including several of the stable 
units (E4-E9). Whereas, for example, femur length (24) 
and tibiotarsus length (26) are next neighbours in the S4 
cluster in BO, in BCM tibiotarus length (26) and tar- 
sometatarsus length (28) are agglomerated at the highest 


level. The larger S1 and S2 units also include some of the 
smaller, topologically stable E units: E4-E6 always oc¬ 
cur within S1, E8-E9 within S3, S4 is always part of S1. 

BOCs identified in the R-cluster analysis of original da¬ 
ta are to a great extent also found in the analysis of body 
mass corrected data (Fig. 3). Again in the body mass cor¬ 
rected data analysis, a main cluster composed of engine¬ 
supporting variables from the flight complex (S2 in BOC 1, 
fig. 3) opposes a main cluster structured by a mixture of 
skeletal elements of the feeding and hind limb locomo¬ 
tion apparatus characters. The main differences to the orig¬ 
inal data analysis are the rather low agglomerative level 
(coefficients of determination in the “size-corrected” 
dataset being two to three times lower than in the origi¬ 
nal dataset), reduction of the involved variable units and 
some changes in variable topology. BOC2, that was com¬ 
posed of SI (including S4), S3, E2, E3 and E7 in the orig¬ 
inal dataset, comprises in the size corrected data set only 
SI (including S4) and E3, the remaining components - 
units S3, E2 and E7 - are clustered at lower agglomera¬ 
tion levels. 

Subcluster SI differs in its hierarchical structure of ag¬ 
glomeration levels of E units and arrangement of variables 
within S clusters. In S2 sternum length (30), keel length 
(32) and keel depth (31) are differently arranged, in S3 
post orbital length (12), and within S4 femur-, tibiotarsus 
and tarsometatarsus lengths (24, 26, 28) changed their po¬ 
sition along branches. 

In both illustrations of single linkage R-cluster analyses 
(Fig. 3), El shows up as a more or less disintegrated unit. 
It consists of measures of ulna (39) and carpometacarpus 
lengths (41), which are tied together at a relatively high 
agglomeration level, especially in the analysis of “size- 
corrected” data, in which both measures are combined 
with carpometacarpus width (42). 

4. DISCUSSION 

4.1. Size and morphological integration 

Do these results confirm or contradict our predictions and 
expectations of hypotheses on morphological integration? 
Corroborating our first assumption, a general body mass 
dependence was identified at a high significance level for 
all individual characters within all three classical function¬ 
al complexes. “Size” could thus be confirmed as a strong 
integrative force in cardueline and fringilline skeletal mor¬ 
phology. Variability was, as assumed in our second hy¬ 
pothesis, relying to a greater extent on body mass in skull 
and beak measurements than in characters allied to the oth¬ 
er functional complexes. Elimination of “size”, decoupled. 
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Table 1 . List of species under study, their body masses and information sources. Sources: DU= Dunning 1993; RE= 
labels of specimens at ZFMK, Bonn; RO= McLean 1988; WA= Ward 2003. 


Species 

body mass in g 

sources 

FringilJa tn on tifringilla 

24 

DU 

Fringilla coelebs 

21,4 

DU 

Coccothraustes coccothraustes 

54 

DU 

Pyrrh it la pyrrh ula 

21,8 

DU 

Pinicola enucleator 

56,4 

DU 

Loxia curvirostra 

40,38 

DU 

Chi oris chloris 

27,8 

DU 

Chloris sinica 

31,3 

DU 

A can this flam mea 

13 

DU 

Carduelis cardnelis 

15,6 

DU 

Linaria cannabina 

15,3 

DU 

Agriospiza jla vi rostris 

15,4 

DU 

Rhodospiza obsoleta 

25,5 

DU 

Spinas spinus 

14,5 

DU 

Spin us cucullatus 

9 

RE 

Spinus barbatus 

16,6 

DU 

Spinus xanthogaster 

12,7 

DU 

Spinus yarrellii 

9,5 

RE 

Spinus psaltria 

9,5 

DU 

Serinus pusillus 

11,6 

DU 

Serinus serinus 

11,2 

DU 

Serinus syriacns 

12,1 

BE 

Serinus canicollis 

13,8 

DU 

Serinus citrinella 

12 

DU 

Alario alario 

11,8 

RO 

Serinus canaria 

12,97 

DU, RE 

Ochrospiza reichenowi 

11,5 

RE 

Ochrospiza atrogularis 

11,4 

DU 

Och rospiza leu copygia 

9 

RE 

Ochrospiza mozambica 

10,6 

DU 

Ochrospiza dorsostriata 

14,4 

DU 

Ochrospiza xanthopygia 

10 

RE 

Crithagra gularis 

15,5 

DU 

Crithagra striolata 

22,4 

DU 

Crithagra albogularis 

25,9 

DU 

Crithagra donaldsoni buchanani 

24 

DU 

Crithagra sulphurata 

19,2 

DU 

Serin ops Jla viven tris 

16,3 

DU 

Crithagra mennelli 

15,3 

RO 

Crithagra burtoni 

29,9 

DU 

Dendrospiza hypostieta 

15 

DU 

Dendrospiza scotops 

15,4 

RO 

Pseudochloroptila totta 

13,4 

WA 
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as predicted in the third assumption, several units of cor¬ 
related variables. Nevertheless, contradicting our expec¬ 
tations, displacement of variables from units or units from 
larger modules was low and robustness of units strong. 
Body mass (size) itself in our sample thus influenced the 
quantity rather than the quality of morphological integra¬ 
tion and did not affect the key feature in fringilline/car- 
dueline skeletal morphology - coevolution of skull and 
hindlimb variables. 

According to classical hypotheses on morphological in¬ 
tegration, “size" is the uppermost integrative level, at deep¬ 
er agglomerative levels, traits covary due to developmen¬ 
tal or functional associations and at still deeper levels due 
to spatial adjacency (Chernoff & Magwene 1999) re¬ 
spectively morphological neighborhood (Alpatov & 
Boschko-Stepanenko 1928). Indeed, in our study, re¬ 
moval of the first hierarchical layer “size", mostly effect¬ 
ed the degree of morphological integration, measurable as 
agglomeration level of character units. Structure and re¬ 
lation of character compositions however remained to a 
great extent constant and robust. For variable units at the 
body mass reduced level, three types of trait correlations 
are identifiable: S-modules are subunits agglomerated at 
higher integrative levels within BOC1 and BOC2, the orig¬ 
inal not body mass reduced, modules. E-modules that are 
not composed of characters from immediate morpholog¬ 
ical neighborhood or partial overlap (E4, E5, E9) can be 
assigned to two categories: Classical functional units (E2, 
E3, E6 and E7) are with one exception (E7) fused at low¬ 
er agglomerative levels (below 0.6 COD) than the units 
of biological roles El and E8 (above 0.8 COD). Our find¬ 
ings thus corroborate the integrative hypotheses quoted 
above, that anatomical association is a primary integra¬ 
tive level in morphology. 

4.2. Interpretation of the observed character mod¬ 
ules 

Eble (2005) classifies four different kinds of modules in 
morphology: 

1. modules due to structural relations, 

2. modules of pleiotropic genotype-phenotype mappings, 

3. developmental units (gene expressions, domains of epi¬ 
genetic dynamics, regions with localized allometric 
growth) and 

4. functional units. 

We would like to extend the definition of the last catego¬ 
ry by modules accomplishing common biological roles 
(Nemesceikal et al. 1992) and thus differentiate between 

4.a. mechanical functioning (like kinetics) and 

4.b. biological functioning (like morphological traits used 

in e. g. feeding). 


Units E4 (postorbital width and intemasal width, 9, 10), 
E5 (skull width and interorbital maximal width. 6, 7) and 
E9 (premaxilla length and total skull length, 11, 16) may 
be explained by their immediate morphological neighbour¬ 
hood and /or partial overlap (“rule of neighbourhood", Al¬ 
patov & Boschko-Stepanenko 1928; structural relations 
Eble 2005). Most of the other results allow for mechan¬ 
ical functional explanations, few for explanations of their 
biological roles. E2 is composed of two skull characters 
(caudal length of the mandible and quadratojugal length, 
3, 13) and a pelvis measure (synsacmm distal length, 23). 
Whereas the first two skull variables are functional coun¬ 
terparts and act in seed husking, their also body mass in¬ 
dependent- correlation with a synsacrum length measure 
comes unawares. E3 combines two mandible depths (4, 
5) with quadratum length (19), characters functioning in 
the kinetism of the avian skull, especially in lateral bill 
movements enhancing seed husking (Nuijens & Zweers 
1997). Also units E6 (skull height and premaxilla depth, 
15, 20) and E7 (femur and tibiotarsus width, 25, 27) are 
easier to explain by constructional constraints than bio¬ 
logical roles. E6 depicts the interdependence of skull and 
beak height, E7 that of leg diameters, in ES.Only units El 
and E8 where measures of the flying apparatus (car- 
pometacarpus (41) and ulna lengths (39)) respectively bill 
tip are tied together (length of pars symphysialis (1) with 
dentary length of the mandible (2)) allow for an interpre¬ 
tation of biological roles: The first unit comprises the 
feather-carrying bones of the wing, responsible for ma- 
noevrability in flight per se. Unit E8 represents “the bill 
tip grasp", a feeding tool well developed in carduelid 
finches (Nemeschkal et al. 1992). 

S4 with the lengths of femur (24), tibiotarsus (26) and tar- 
sometatarsus (28) obviously reflects a character pattern of 
concerting limb dimensions as guided by developmental 
processes during ontogeny (Nemeschkal 1999). 

Comparison of interspecific correlation patterns 
(macroevolutionary pattern, this study Fig. 3) with infra- 
specific analyses of variable correlations (microevolution¬ 
ary level, Nemeschkal & van den Elzen 1994, fig. 7) re¬ 
veal that six modules (E3, E4, E5, E7, E8 and E9 of this 
study) correspond to units also found within carduelid 
finches and three modules (E5, E7 and E9) even to units 
occurring within the phylogenetically distant pigeons. At 
the microevolutionary level Nemeschkal et al. (1992) pro¬ 
vided evidence for clade specificity of variable units (mod¬ 
ules) and their correspondence to expressions of develop¬ 
mental control genes (Nemeschkal 1999). Morphologi¬ 
cal modularity might generally be seen as product of de¬ 
velopmental modularity, because morphological patterns 
of organization emerge in ontogeny (Eble 2005). Conse¬ 
quently, each of the extracted variable complexes might 
additionally also correspond to developmental or body 
plan modules. 
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Table 2. Skeletal variables under study and their coefficients of determination. BO represent original data, BCM resi¬ 
duals for the body mass reduced variables. 


Skeletal variables 


Coeff. of determination 


Skull measures 

01 mandible, pars symphysialis length 0,7014 

02 mandible, dentary length 0,7641 

03 mandible, pars caudalis length 0,7642 

04 mandible, dentary depth 0,6924 

05 mandible, pars caudalis depth 0,7883 

06 skull width 0,8239 

07 interorbital, maximal width 0,8353 

08 interorbital, minimal width 0,6704 

09 postorbital width 0,8048 

10 intemasal width 0,7856 

11 premaxilla length 0.8240 

12 postorbital length 0,86 

13 quadratojugal length 0,7953 

14 orbital width 0,7888 

15 skull height 0,7280 

16 skull length 0,8610 

17 palatinum width 0,6176 

18 pterygoid length 0,6932 

19 processus orbitalis quadrati length 0,7988 

20 premaxilla depth 0,7911 


Pelvic girdle and hindlimbs 

21 synsacrum length 

22 synsacrum width 

23 synsacrum distal length 

24 femur length 

25 femur width 

26 tibiotarsus length 

27 tibiotarsus width 

28 tarsometatarsus length 

29 tarsometatarsus width 


Pectoral girdle and forelimbs 

30 sternum length 0,8413 

31 keel depth 0,6897 

32 keel length 0,7632 

33 sternum, distal width 0,6617 

34 sternum depth 0,7393 

35 coracoid length 0,8401 

36 scapula length 0,8114 

37 humerus length 0,7644 

38 humerus, proximal end width 0,6246 

39 ulna length 0,3233 

40 ulna width 0,6235 

41 carpometacarpus length 0,2333 

42 carpometacarpus width 0,6454 


0,7748 

0,7076 

0,7464 

0,8355 

0,7781 

0,7675 

0,7562 

0,6180 

0,3034 
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How genes and developmental pathways co-operate is best 
documented for the avian beak. It is constructed by mul¬ 
tiple facial prominences (Helms et al. 2005, Helms & 
Brugmann 2007), where the frontal nasal mass (FNM), 
lateral nasal prominences (LNP), and maxillary promi¬ 
nences (MXP) comprise the upper beak and the mandibu¬ 
lar prominence (MDP) forms the lower beak. The identi¬ 
ties of facial prominences are specified early in the neu¬ 
ral crest stage and may involve homeobox genes such as 
Hox (Couly et al. 1998) or MSX (Brown et al. 1997, Wu 
et al. 2006). 

In metazoa, at least 17 signal transduction pathways op¬ 
erate to activate or repress different genes at distinct times 
and places in the embryo. Five predominate in early em¬ 
bryonic development: the Wnt, TGF-beta, Hedgehog, 
RTK, and Notch pathways. Five more are used in late de¬ 
velopment, and seven more in the functions of differen¬ 
tiated cells (Gerhart 1999). 

Diversity in beak shapes of Galapagos finches, e. g., is 
known to originate from the activity of a morphogenetic 
bone protein (BMP4 of the TGF-beta family; Abzhanov 
et al. 2004) and calmodulin, a protein that binds and trans¬ 
ports calcium ions (CaM; Abzhanov et al. 2006). Where¬ 
as BMP4 stimulated growth of beaks along two dimen¬ 
sions, - it produced deeper and wider beaks and thus ex¬ 
plained the linkage in the variation of these traits, - the 
authors found CaM as regulator of beak length. Thus in 
Galapagos finches, beak length develops independently 
from width and depth due to two different factors that lead 
to changes in growth along different dimensions 
(Abzhanov et al. 2006). In ducks, chickens and cockatiels, 
BMP4 is also differently expressed (Wu et al. 2004, Wu 
et al. 2006). At late stages of development, chicken and 
duck embryos had two localized growth zones in the 
FNM, which melted in the chicken, but stayed separated 
in the duck. Ducks, moreover exhibited a wider FNM and 
more activity in another growth factor (fibroblast growth 
factor 8). In cockatiels a thicker FNM increased in a dif¬ 
ferent direction and the mandibular prominence (MDP) 
was suppressed. BMP4 was involved in all species in me¬ 
diating activity in all localized growth zones. Experimen¬ 
tal overexpression in BMP4 altered beak shapes among 
all species and beak curvation was induced by asymmet¬ 
ric growth activity in a facial prominence. BMP4 was al¬ 
so found to be responsible for regulating a homeobox 
gene: As it increased the expression of MSX 1, the authors 
suggested that MSX1 activity is regulated by the BMP 
pathway. 

4.3. Conclusions 

Under the assumption that BMP4 is a main regulator for 
the expression of developmental control genes in the avian 


skeleton, this protein, for instance, can be seen as a key 
growth factor accounting for size and shape variation in 
our skeletal variables. A developmental growth factor can 
also explain the strong size dependence of single variables, 
but cannot elucidate the interspecific correlation patterns 
of variable units. Thus in the light of current knowledge 
on gene expression variability and research in develop¬ 
mental genetics the classical view on morphological in¬ 
tegration and genesis of modularity might be too one¬ 
sided. 

Wagner et al. (2007) offer a broader approach to modu¬ 
larity. Summarizing empirical evidence, they discern be¬ 
tween two levels: The elements the modules consist of and 
their kinds of interaction. Elements vary from molecular 
level (nucleotides) to morphological traits and their con¬ 
nections from physical to dynamical and statistical. They 
define three kinds of modules: Variational, functional and 
developmental modules, but stress a duality, as modules 
may belong to several categories at the same time. Units 
observed in the present study on cardueline/fringilline 
skeletal variables fall in the category of variational mod¬ 
ules, as their measure is statistical. For the interpretation 
of variational modules the authors offer several answers: 
identical developmental origin, similarity of covariation 
and gene expression territories, and pleiotropic effects. So 
our variable units (or character complexes; Nemeschkal 
et al. 1992) can be defined not only as a duality but a plu¬ 
rality: according to their operability they are functional 
modules, according to their mode of exploration they are 
variational modules, and according to their genesis devel¬ 
opmental modules. 

The open problem is, whether modules arise through the 
action of natural selection or because of biased mutation¬ 
al mechanisms (Wagner 2007). Both neutral models and 
models based on natural selection are offered: e. g. Lynch 
(2007) favors the neutral model and stated that “... emer¬ 
gent biological features such as complexity, modularity, 
and evolvability, all of which are current targets of con¬ 
siderable speculation, may be nothing more than indirect 
by-products of processes operating at lower levels of or¬ 
ganization.’ 1 We do not follow the author in refusing nat¬ 
ural selection and external evolutionary forces, but as¬ 
sume, that selection via fitness promotes certain pheno¬ 
types derived from their genetic basis and processes in on¬ 
togeny. Under the viewpoint of developmental constraints 
as internal evolutionary forces, we like to generalize the 
findings of Wu et al. (2004, 2006) on the origin of beak 
shape in birds. The authors offer a more reconciling con¬ 
clusion, that morphological diversity may be achieved by 
modulating prototypical molecular and cellular modules. 
Transforming growth factors may mediate the range, lev¬ 
el, and duration of locally enhanced growth, thus provid¬ 
ing a spectrum of morphological designs for selection (Wu 
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et al. 2006). Also Wagner et al. (2007) arrive at the con¬ 
clusion that mutational processes (internal processes) fa¬ 
vor the origin of modularity and selection pressures (ex¬ 
ternal processes) reinforce the mutational bias. 
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